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The origin and nature of the highest energy cosmic ray events is currently the subject of intense investigation 
by giant air shower arrays and fluorescent detectors. These particles reach energies well beyond what can be 
achieved in ground-based particle accelerators and hence they are fundamental probes for particle physics as well 
as astrophysics. One of the main topics today focuses on the high energy end of the spectrum and the potential 
for the production of high-energy neutrinos. Above about 10 20 eV cosmic rays from extragalactic sources are 
expected to be severely attenuated by pion photoproduction interactions with photons of the cosmic microwave 
background. Investigating the shape of the cosmic ray spectrum near this predicted cut-off will be very important. 
In addition, a significant high-energy neutrino background is naturally expected as part of the pion decay chain 
which also contains much information. 

Because of the scarcity of these high-energy particles, larger and larger ground-based detectors have been built. 
The new generation of digital radio telescopes may play an important role in this, if properly designed. Radio 
detection of cosmic ray showers has a long history but was abandoned in the 1970's. Recent experimental develop- 
ments together with sophisticated air shower simulations incorporating radio emission give a clearer understanding 
of the relationship between the air shower parameters and the radio signal, and have led to resurgence in its use. 
Observations of air showers by the SKA could, because of its large collecting area, contribute significantly to 
measuring the cosmic ray spectrum at the highest energies. Because of the large surface area of the moon, and 
the expected excellent angular resolution of the SKA, using the SKA to detect radio Cherenkov emission from 
neutrino-induced cascades in lunar regolith will be potentially the most important technique for investigating 
cosmic ray origin at energies above the photoproduction cut-off. 



1. Observational and theoretical motiva- 
tion 

Understanding the origin of the ultra-high en- 
ergy cosmic rays (UHECR), the highest energy 
particles observed in nature, is of great impor- 
tance as it may impact our understanding of 
particle physics, fundamental cosmology, and ex- 
tremely energetic phenomena in the Universe. 
The energy spectrum of UHECR extends up to at 
least 10 11 GeV, and in the rest frame of a UHECR 
proton, photons of the 2.73 K cosmic microwave 
background radiation (CMBR) are strongly blue- 
shifted to gamma-ray energies. The threshold for 
Bcthe-Heitler pair production and pion photopro- 



duction by UHECR protons on the CMBR are 
close to 2 x 10 8 GeV and 2 x 10 10 GeV, such that 
protons at 3xl0 10 GeV and 3xlO n GeV typi- 
cally lose a large fraction of their energy in a time 
of 1 Gpc/c (3 x 10 9 y) and 10 Mpc/c (3 x 10 7 y), 
respectively. This would imply that sources of 
ultra-high energy cosmic rays would have to be 
close if the particles themselves behave as pre- 
dicted. The importance of pion photoproduction 
on the CMBR was first noted by Greisen (1966) 
and Zatsepin & Kuzmin (1966) and the cut-off 
they predicted at ~ 10 11 GeV is referred to as 
the "GZK cut-off" . 

Yet, until today neither the nature of the par- 
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tides nor of their accelerators has been revealed. 
It is well established that in some astrophysical 
magnetized plasma regions particles (leptons) are 
accelerated, but whether and where this holds for 
UHECR is unclear. Some basic constraints can 
nonetheless be given. Cosmic ray acceleration 
sites must be large enough to contain the gyrora- 
dius of the accelerated particles, as well as hav- 
ing scattering centres with appropriate velocities. 
In addition, the acceleration must be sufficiently 
rapid that high energies can be achieved in an 
accelerator's lifetime, and that energy losses by 
pion photoproduction and synchrotron radiation 
do not cut off the spectrum too soon (Hillas 1984). 
It is currently unknown whether the UHECR 
are Galactic or extragalactic in origin. Compo- 
sition measurements are also important because 
if UHECR are observed to include nuclei other 
than protons then these must be from Galactic, 
or very nearby extragalactic acceleration sources 
to avoid photodisintegration (see e.g. Yamamoto 
et al. 2004). However, the promising extragalac- 
tic source candidates for UHECR above 10 11 GeV 
are typically at distances too far for UHECR 
to reach us unaffected by interactions with the 
CMBR. This is the basic dilemma we are faced 
with today. 

1.1. UHECR observations 

Below the GZK cut-off UHECR may, to some 
extent, point back to their sources depending 
on the structure and strength of the magnetic 
field between the sources and our Galaxy. No 
statistically compelling anisotropy has been de- 
tected in the UHE CR. The energy spectrum of 
UHECR detected by AGASA is shown in Fig. □ 
There are two main problems at present: the 
flux of UHECR is so low that few events have 
been detected for reliable conclusions concerning 
the presence or absence of a GZK cut-off or any 
anisotropy, and in the case of a GZK cut-off spec- 
tral information above the cut-off would be lost. 
New experiments such as HiRes and the 3000 km 2 
Peirre Auger Observatory will help to address the 
the first issue and, because of its huge area, use 
of the SKA may also help here by direct radio de- 
tection of UHECR air showers. The second prob- 
lem, that of loss of spectral information above the 
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Figure 1. Spectrum of UHECR detected by 
AGASA. Numbers attached to points show the 
number of events in each energy bin. (From 
Takeda et al. 2003.) 



GZK cut-off, is best explored through UHE neu- 
trino astronomy and again the use of the SKA, 
either to detect neutrino-induced air showers or 
radio Cherenkov bursts from electromagnetic cas- 
cades in lunar regolith initiated by interactions of 
UHE neutrinos, has the potential to greatly add 
to our knowledge of the the origin of the highest 
particles in nature. 

1.2. The GZK problem 

Due to interactions with the CMBR, there is 
expected to be a spectral downturn, the GZK 
cut-off, for particles which have travelled more 
than a few tens of Mpc. Figure |21 shows the 
distribution in energy as a function of distance 
travelled for UHECR protons with initial ener- 
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Figure 2. (a) Arrival energy distribution for 
protons injected with energy between 10 21 ' 9 and 
10 22 eV after propagation on 10, 20, ... 200 
Mpc. (b) Fractional energy contained in nucleons 
(solid line), 7-rays from photoproduction (long 
dashes) and BH pair production (short dashes) 
for protons injected with a E~ 2 power law spec- 
trum with an exponential cutoff at 10 215 eV. The 
dash-dot lines show the fractional energy in muon 
(long) and electron (short) neutrinos and antineu- 
trinos. (From Stanev at al. 2000.) 



gies close to 10 22 eV. As can be seen, after 80 Mpc 
no protons have energies above 3 x 10 20 eV. Of 
course, UHECR protons are also deflected by ex- 
tragalactic magnetic fields, and so any source of 
3 x 10 20 eV UHECR would need to be much nearer 
than 80 Mpc. However, several experiments have 
reported CR events with energies above 10 20 eV 
with the highest energy event having 3 x 10 20 eV 



(Bird et al.1995). Very recent data from the two 
largest aperture high energy cosmic ray detectors 
are contradictory: AGASA (Takeda et al. 2003) 
observes no GZK cut-off while HiRes (Abbasi et 
al. 2004) observes a cut-off consistent with that 
expected. A systematic over-estimation of energy 
of about 25% by AGASA or under-estimation of 
energy of about 25% by HiRes could account the 
discrepancy (Abbasi et al. 2004), but the con- 
tinuation of the UHECR spectrum to energies 
well above 10 20 eV is now far from certain. Fu- 
ture measurements with Auger (Auger Collabora- 
tion 2001) should resolve this question. Whether 
or not the spectrum does extend well beyond 
10 20 eV, determining the origin of these parti- 
cles could have important implications for astro- 
physics, cosmology and particle physics. 

1.3. The acceleration problem 

By plotting magnetic field vs. size of various 
astrophysical objects (Fig.|3Jl, Hillas (1984) iden- 
tified possible sites of acceleration of UHECR 
based on whether the putative source could con- 
tain the gyroradius of the accelerated particles, 
and on the likely velocity of scattering centres 
in these sites. Following Hillas (1984) one finds 
that possible sites included neutron stars (10 7 - 
10 13 G), gamma ray bursts and active galactic nu- 
clei (10 3 -10 4 G), and lobes of giant radio galaxies 
and galaxy clusters (10 _7 -10 _5 G). This identifi- 
cation of possible sources does not take account 
of energy losses (synchrotron) and interactions 
(Bethe-Heitler and pion photoproduction) which 
can cut off the spectrum, and so apply an addi- 
tional constraint which we discuss below (in his 
original paper Hillas (1984) also used this addi- 
tional constraint to narrow the field of possible 
sources). 

For particle acceleration by electric fields in- 
duced by the motion of magnetic fields B (in- 
cluding those at astrophysical shocks), the maxi- 
mum rate of momentum gain by relativistic par- 
ticles of charge Ze can be written (in SI units) 
(dp/dt) acc = £(p)ZecB where £(p) < 1 is the ac- 
celeration rate parameter and depends on the de- 
tails of the acceleration mechanism (see the re- 
view by Jones & Ellison 1991, on the plasma 
physics of shock acceleration, which also includes 



4 




Figure 3. "Hillas plot" showing (chain curves) 
magnetic field vs. gyroradius for proton momenta 
10 15 , 10 16 , 10 24 eV/c. The solid curves 

bound the parameter space of accelerated parti- 
cles for a given acceleration rate parameter (see 
text). Typical size and magnetic field of pos- 
sible acceleration sites (taken from Hillas 1984) 
are shown for neutron stars (ns), white dwarfs 
(wd), sunspots (ss), magnetic stars (ms), AGN 
(ag), interstellar space (is), supernova remnants 
(sn), radio galaxy lobes (rg), galactic disk (d) 
and halo (h), clusters of galaxies (cl) and in- 
tergalactic medium (ig). Typical jet-frame pa- 
rameters of the synchrotron proton blazar model 
(Miiecke et al. 2003) and gamma ray burst model 
(Pcllctier & Kersate 2000) are indicated by the 
open squares labelled "bl" and and "gb" . (From 
Protheroe 2004.) 



a brief historical review and refers to early work) . 
To estimate cut-off momenta (or energies), one 
needs the acceleration rate. The following val- 
ues for the acceleration rate parameter have been 
suggested: maximum possible acceleration rate 
£(Pcut) = l, plausible acceleration at perpendicular 
shock with speed 0.1c, £(f> C ut) ~ 0.04, and plausi- 
ble acceleration at parallel shock with speed 0.1c, 



Figure 4. Maximum energy as a function of mag- 
netic field of protons for maximum possible accel- 
eration rate £ = 1 (upper solid curve), £ = 0.04 
(middle solid curve), £ = 1.5 x 10~ 4 (lower solid 
curve). Dashed curves are limits from Bcthc- 
Heitler pair production and pion photoproduction 
only (solid curves include synchrotron loss). Dot- 
dot-dot-dash curves are lines of constant Larmor 
radius as labelled. (From Protheroe 2004.) 



£(Pcut) ~ 1.5 xlO -4 (Protheroe 2000). Based on 
the total momentum loss rate for Bethe-Heitler 
pair production and pion photoproduction on the 
CMBR, synchrotron losses and redshifting the 
proton cut-off momentum is plotted in Fig. 0] as 
a function of magnetic field for three adopted £- 
values (chain lines arc for constant Larmor ra- 
dius as labelled). This plot clearly shows that 
to accelerate protons to ~ 10 20 eV large regions 
of relatively low magnetic field ~ 10 _7 -10 _3 G 
are needed, apparently ruling out high magnetic 
field regions for the origin of UHECR (see also 
Medvedev 2003). One sees that, in principle, pro- 
tons can be accelerated up to ~ 5 x 10 22 eV in 
Mpc scale region with ~ 10" 5 G. 

Returning to the Hillas plot (Fig. [3]), con- 
straints have been added corresponding to the 
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three curves in Fig.^J and the chain lines give con- 
stant proton energy values as indicated. Sources 
to the right of the solid curves are excluded; a pos- 
sible exception to this is in the case of relativisti- 
cally beamed sources (e.g. for AGN see Protheroe 
et al. 2003, and for GRB see Pelletier & Kersate 
2000) where neutrons emitted along the direction 
of relativistic motion can be Doppler boosted sig- 
nificantly in energy. Another possible exception 
is the case of so called "one-shot" mechanisms 
(e.g. Haswell et al . 1992, Sorrcll 1987) where a 
particle is accelerated by an electric field along 
a nearly straight path which is essentially paral- 
lel to the magnetic field such that curvature and 
synchrotron losses are negligible. Suggested sites 
for this include polarization electric fields arising 
in plasmoids injected into a neutron star's mag- 
netosphere (Litwin & Rosner 2001) and magnetic 
re-connection in the magnetosphere of accretion 
induced collapse pulsars (de Gouvela Dal Pino 
& Lazarian 2001). Another possibility is plasma 
wakefield acceleration, i.e. acceleration by collec- 
tive plasma waves, possibly in the atmosphere of a 
GRB, or "surf-riding" in the approximately force- 
free fields of the relativistic wind of a newly born 
magnetar (Arons 2003). In these cases it is un- 
clear whether the requirements of negligible radi- 
ation losses can be met. 

Alternative scenarios for UHECR origin in- 
clude emission and decay of massive particles 
( "X-particles" ) by topological defects (TD) or de- 
cay of massive primordial particles. Because of 
the resulting flat spectrum of particles (including 
neutrinos, gamma-rays and protons) extending 
possibly up to GUT (grand unified theory) scale 
energies, topological defect models have been in- 
voked to try to explain the UHECR. Propagation 
of the spectra of all particle species over cosmolog- 
ical distances is necessary to predict the cosmic 
ray and gamma-ray spectra expected at Earth. 
In most cases this results in excessive gamma- 
ray fluxes at GeV energies in addition to cosmic 
rays. Massive relic particles on the other hand, 
would cluster in galaxy halos, including that of 
our Galaxy, and may give rise to anisotropic cos- 
mic ray signals at ultra high energies. One pos- 
sibility proposed for getting around this is if an 
extragalactic source emits a very high luminos- 



ity in UHE neutrinos, some of which interact 
with relic neutrinos gravitationally bound to our 
galaxy producing "Z-bursts" which generate the 
events observed above the expected GZK cut-off. 
(see Protheroe & Clay 2004 for a recent review of 
UHECR.) 

1.4. Neutrino signatures 

One way of getting information about accel- 
eration sources of UHECR is through the spec- 
tral shape near acceleration cut-off. One of the 
present authors (Protheroe 2004) has recently 
shown that in the case of protons the spectrum 
can actually be quite sensitive to the astrophys- 
ical acceleration environment. Despite the fact 
that for extragalactic UHECR almost all spec- 
tral information above the GZK cut-off is lost, 
significant information is preserved in the spec- 
trum of neutrinos produced as a result of pion 
photoproduction interactions during propagation 
(Protheroe 2004). Furthermore, the spectrum 
of these GZK neutrinos differs significantly from 
that in Z-burst and topological defect (TD) sce- 
narios, and of course the neutrinos are not de- 
flected by magnetic fields and so should point 
back to where they were produced. Hence, UHE 
neutrino astronomy will be able provide much 
needed clues to the origin of the UHECR. 

The spectra of protons and neutrinos escaping 
from an acceleration source and after propagation 
for 100 Mpc/c is shown in Fig.0for various possi- 
ble acceleration environments represented by the 
spectrum of magnetic turbulence present (power- 
law dependence of acceleration rate), the average 
magnetic field, its alignment, speed of scattering 
centres (acceleration rate and maximum energy) , 
and the size of the acceleration region (decay or 
escape of photoproduced neutrons). As can be 
seen, while there is little difference in the spec- 
trum of UHECR after propagation over 100 Mpc, 
much information is preserved in the spectrum of 
UHE neutrinos ( "GZK neutrinos" ) produced dur- 
ing propagation as the UHECR flux is eroded by 
the GZK-cutoff effect. Of course, for very dis- 
tant sources UHECR would not be expected to 
be observed from directions of sources. Indeed, 
few if any may arrive at all because of difficulty 
in reaching Earth through extragalactic magnetic 
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Figure 5. Spectra of protons and neutrinos (all flavours) escaping from the acceleration region (dotted 
curves) and after propagation for time 100 Mpc/c (solid curves) forp cut c = 10 20 (leftmost curves), 10 20 5 , 
. . . , 10 23 eV (rightmost curves); T = 2. (From Protheroe 2004.) 



fields, whereas UHE neutrinos will arrive essen- 
tially undeflected. A very sensitive UHE neu- 
trino telescope may therefore observe neutrinos 
from extragalactic UHECR sources. The diffuse 
GZK neutrino background can actually be quite 
large if the UHECR sources evolve strongly with 
rcdshift (Engel ct al. 2001). Nevertheless, huge 



collecting areas will be required for the detection 
of UHE neutrinos, and it is here that the SKA 
through direct detection of neutrino-induced air 
showers and, perhaps more importantly, through 
the detection of Chcrcnkov radio transients from 
neutrino-induced showers in lunar regolith may 
make a major contribution to understanding the 
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origin of the UHECR. 

2. Detection of high energy particles: his- 
torical background 

There is clearly strong theoretical and phe- 
nomcnological motivation to detect both the pre- 
sumably hadronic cosmic rays and associated 
neutrinos at EeV to ZeV energies. The difficulty 
arises from the extremely low fluxes present-for 
the highest energy cosmic rays at or above the ~ 
6 x 10 19 eV GZK cutoff, one can expect of order a 
few per km 2 per century at most. The associated 
EeV neutrino fluxes are, in the most optimistic 
scenarios, perhaps 1-2 orders of magnitude larger 
than this, but their detection efficiency is at most 
~ 1% per cubic km of water- equivalent material, 
and thus the neutrino rates are abysmally low in 
all existing and most planned detectors (though 
a recently approved NASA long-duration balloon 
experiment, ANITA (Barwick et al. 2003), may 
get an early, low-resolution view of these fluxes). 

2.1. Giant air shower detectors. 

Since the early 1960's through the mid-1980's 
the highest energy cosmic ray detectors were ex- 
clusively large ground arrays of scintillators or 
Cherenkov counters making direct detection of 
secondary particles, mainly electrons, gamma- 
rays, and muons within the confines of the air 
shower itself as it impacts the ground (see Nagano 
and Watson 2000 for a review and references to 
major air shower detectors). At the highest ener- 
gies, air shower detectors gain much of their col- 
lecting aperture by capturing the edges of show- 
ers whose cores fall outside their fiducial array 
boundaries, sometimes by many hundreds of me- 
ters. Thus the shower energy must be estimated 
by parametric models for the particle density at 
the shower periphery - a technique which has un- 
dergone much evolution throughout the history of 
giant air shower detection, and still retains much 
controversy in the details of its application even 
today. 

In the mid-1980's the first Nitrogen air fluo- 
rescence (N2fl) detector, the Fly's Eye, came on- 
line. Since that time, both the Fly's Eye and the 
follow-on High Resolution Fly's Eye (HiRes) have 



become competitive with the air shower ground 
arrays in their detection efficiency and aperture, 
and HiRes now has the largest exposure and data 
sample of any detector to date. 

The N2fl technique is very different than that 
of ground array detection, since the detectors 
do not require direct intersection with any por- 
tion of the air shower particles, but rather de- 
tect the secondary incoherent radiation from de- 
excitation of Nitrogen heated by the passage of 
the shower. Such emission may be seen by opti- 
cal telescopes of several m 2 aperture out to tens 
of km distance from the shower itself, and thus a 
small installation of modest, low-optical quality 
(e.g., searchlight-style) reflectors can, by viewing 
a good fraction of the surrounding sky, create an 
effective air shower collecting aperture of several 
thousand km 2 sr. The only drawbacks to this 
technique are its sensitivity to atmospheric atten- 
uation in the near ultraviolet (where the nitrogen 
emission lies), and its requirement for complete 
darkness and clear weather. These constraints 
lead to a low net long-term duty cycle of less 
than 10%, compared to the ~ 100% duty cycle 
of a ground array. 

2.2. Cosmic ray air shower radio detection 

Interest in radio techniques for giant air shower 
detection stemmed originally from the suggestion 
by Askaryan (1962) that any electromagnetic cas- 
cade in a dielectric material (gas, liquid or solid) 
should rapidly develop net negative charge asym- 
metry due to electron scattering processes and 
positron annihilation. The net electronic charge 
excess was estimated to be ~ 20 — 30%, and 
Askaryan proposed that Cherenkov radiation at 
wavelengths larger than the longitudinal dimen- 
sions of the shower (~ 1 m in air, and ~ 1 cm 
in liquids or solids) would be emitted coherently, 
yielding a quadratic scaling of received power 
with the shower energy. This latter property 
immediately suggests that radio emission might 
dominate the secondary radiation at the high- 
est energies. We defer discussion of this so-called 
Askaryan effect in solids to a later section; how- 
ever, its application to air showers was immedi- 
ately noticed and pursued. 
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2.2.1. History 

Radio emission from cosmic ray air showers was 
discovered for the first time by Jelley and co- 
workers in 1965 at a frequency of 44 MHz. They 
used an array of dipole antennas in coincidence 
with Geiger counters. The results were soon veri- 
fied and emission from 2 MHz up to 520 MHz was 
found in a flurry of activities in the late 1960's. 
These activities ceased almost completely in the 
subsequent years due to several reasons: difficulty 
with radio interference, uncertainty about the in- 
terpretation of experimental results, and the suc- 
cess of other techniques for air shower measure- 
ments. 

The radio properties of air showers are sum- 
marized in an excellent and extensive review by 
Allan (1971). The main result of this review can 
be summarized by an approximate formula relat- 
ing the received voltage of air showers to various 
parameters, where we also include the presumed 
frequency scaling: 



e„ = 20 /iV m" 1 MHz" 1 ( 1Q ^ P y I *»> " 



cos 9 exp 
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(1) 



Here E p is the primary particle energy, R is the 
offset from the shower center and Rq is around 
110 m at 55 MHz, is the zenith angle, a is the 
angle of the shower axis with respect to the geo- 
magnetic field, and v is the observing frequency 
(see also Allan et al. 1970; Hough & Prescott 
1970). The leading factor of 20 has been disputed 
over the years since it was first published, and 
could be an order of magnitude smaller. 1 

The voltage of the unresolved pulse in the co- 
herent regime (y < 100 MHz) can be converted 
into an equivalent flux density (a flux density for 
a steady continuum source required to produce 
the same energy over the bandwidth limited time 
interval At) in commonly used radio astronomical 
units 



1 More likely, the controversy over this coefficient probably 
stems from the wide variation in measurement conditions 
and the uncertainties in the flux calibration of the radio 
antennas as well as in the energy calibration of the parti- 
cles. 
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The pulse duration is At ~ 1/ Av if the mea- 
surement is bandwidth-limited. Note, that for 
larger bandwidths and hence higher time resolu- 
tion the energy of the pulse itself does not in- 
crease, however, the equivalent flux density of a 
steady source needs to increase, in order to pro- 
duce an energy comparable to the pulse in the 
shorter time interval. In the earlier measurements 
the pulses were always unresolved when observing 
with Av ~ 1 MHz. 

The formula was determined experimentally 
from data in the energy regime 10 16 eV < 
E p < 10 18 eV. The flux density around 100 MHz 
seems to depend on primary particle energy as 
S v oc E'l (Hough & Prescott 1970; Vernov et al. 
1968; Fig. EJ as expected for coherent emission 
(see below). This dependency is, however, not 
yet undoubtedly established, since a few earlier 
measurements apparently found somewhat flat- 
ter power-laws (Barker et al. 1967 as quoted in 
Allan 1971). 

Very little concrete data exist on the spectral 
dependence of EAS (Extensive Air Shower) radio 
emission (e.g., Spencer 1969). Figure shows a 
tentative EAS radio spectrum with a v~ 2 depen- 
dence for the flux density (y~ x dependence for the 
voltage). The 2 MHz data point was made with a 
different experiment and there is a possibility that 
the spectrum is somewhat flatter between 10-100 
MHz (see Datta et al. 2000), but this is not ver- 
ified. The polarisation of the emission could be 
fairly high and is basically along the geomagnetic 
E-W direction (Allan, Neat, & Jones 1967) which 
strongly supports an emission mechanism related 
to the geomagnetic field. 

Finally, one needs to consider the spatial struc- 
ture of the radio pulse. The current data 
strongly supports the idea that the emission is 
not isotropic but is highly beamed in the shower 
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Figure 6. The dependence of EAS radio flux on 
the primary particle energy as measured by Ver- 
nov et al. (1968) following roughly a power- 
law. Some earlier papers found somewhat flatter 
dependencies. 



direction. Figure [8] shows EAS radio pulse am- 
plitude measurements function of distance 
R from the shower axis - the flux density drops 
quickly with offset from the center of the shower. 
The characteristic radius of the beam is of order 
100 meter for a 10 17 eV vertical shower, with the 
emission originating at 5-7 km distance above an 
observer at sea level. The implied angular diam- 
eter of the beam is thus 9 ~ 0.2/6 = 1.9°. 

2.2.2. The synchrotron model and recent 
work. 

Experiments have clearly established that cos- 
mic ray air showers produce radio pulses. The 
original motivation was due to a suggestion from 
Askaryan (1962) who argued that annihilation of 
positrons would lead to a negative charge excess 
in the shower, thus producing Cherenkov radi- 
ation as it rushes through the atmosphere. At 
radio frequencies the wavelength of the emission 
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Figure 7. A tentative radio pulse spectrum for 2 
MHz to 520 MHz. The data are not simultaneous. 
From Allan (1971) and Spencer (1969). 



is larger than the size of the emitting region and 
the emission should be coherent. The radio flux 
would then grow quadratically with the number 
of particles rather than linearly and thus would 
be greatly enhanced. This effect is important in 
dense media where it was already experimentally 
verified (Saltzberg et al. 2001; see below) and is 
important for detecting radio emission from neu- 
trino showers in ice or on the moon. 

However, the dependence of the emission on 
the geomagnetic field detected in several later ex- 
periments indicates that another process may be 
important. The basic view in the late 60's was 
that the continuously created electron-positron 
pairs were then separated by the Lorentz force 
in the geomagnetic field which led to a trans- 
verse current in the shower. If one considers a 
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Figure 8. Normalized radio pulse amplitudes in 
/LtV m _1 MHz -1 at 55 MHz as a function of dis- 
tance R in meters from the shower axis. Each 
data point corresponds to one measured cosmic 
ray event. The amplitudes were normalized to a 
reference energy of E p = 10 17 eV assuming the 
above mentioned linear dependence of voltage on 
primary particle energy. The measurements were 
made for zenith angles 9 < 30°. Crosses and dots 
represent different particle energy bins between 
10 17 eV and 10 18 eV. The plus sign at 500 meters 
marks a single 10 19 eV event. From Allan (1971). 



frame moving along with the shower, one would 
observe electrons and positrons drifting in oppo- 
site directions impelled by the transverse electric 
held induced by the changing geomagnetic flux 
swept out by the shower front. (Only in the case 
of shower velocity aligned with the magnetic held 
lines will this induced electric held vanish). This 
transverse current then produces dipole (or Lar- 



mor) radiation in the frame of the shower. When 
such radiation is Lorentz-transformed to the lab 
frame, the boost then produces strongly forward- 
beamed radiation, compressed in time into an 
electro-magnetic pulse (EMP). This was calcu- 
lated by Kahn & Lerche (1966) and also Colgate 
(1967). 

Falcke & Gorham (2003) suggested it might 
be better to think of the emission simply as be- 
ing synchrotron-like in the earth's magnetic held, 
or "coherent geosynchrotron emission" , as they 
called it. This process is probably equivalent to 
the previous suggestions since it is derived from 
the basic formula for dipole radiation and the 
Poynting vector but does not require a consider- 
ation of charge separation: The different sign of 
the charges is canceled by the opposite sign in the 
Lorentz force for electrons and pairs and hence 
both contribute in exactly the same way to the 
total flux (radio astronomers will surely remem- 
ber that an electron/positron plasma produces al- 
most the same amount of synchrotron emission as 
a pure electron plus proton plasma). 

The basic and intuitive derivation of this ef- 
fect can be found in Falcke & Gorham (2003) us- 
ing standard synchrotron radiation theory. One 
important effect which is explicitly neglected by 
this simple treatment is the Fresnel zone prob- 
lem - vertical air showers at 10 19 eV reach their 
particle maximum at ground level, and the radio 
emission arrive nearly simultaneously to the par- 
ticle "pancake," indicating that the far-field con- 
ditions, where the radiation field has had time 
to become well-separated from its source, are not 
satisfied. Any estimate of the details of the re- 
ceived radio emission which is intended to help 
with detailed detector design, such as what may 
be required to justify any impact on SKA parame- 
ters or planning, must therefore treat the problem 
with much greater fidelity. 

Such high-fidelity simulations of geosyn- 
chrotron emission are now beginning to appear 
in the literature, and as the interest in this ap- 
proach grows, along with the compelling nature 
of the ultra-high energy cosmic ray problem, the 
simulations can be expected to improve as well. 
In the following section, we describe recent results 
in this direction. 
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2.2.3. Air shower electrodynamics: de- 
tailed modeling and Monte Carlo 
simulation. 

The challenge of developing high-fidelity air 
shower radio simulations breaks into three dis- 
tinct problems: 

1. The adaptation of existing air shower sim- 
ulation codes to provide the particle identi- 
fication and sampling needed for electrody- 
namics modeling; 

2. The implementation of actual electrody- 
namics computation within the modified air 
shower code, and the development of radi- 
ation propagation model; and 

3. the modeling of the detector geometry and 
detection process. 

To date, no group has implemented all three 
aspects of this program, but we describe here two 
efforts which have gone much further than others 
in addressing the difficult problem of the electro- 
dynamics and detection modeling. 

Simulations by the Chicago/Hawaii 
group. 

One result with a first-order electrodynamics 
Monte Carlo simulation has been completed by 
Suprun et al. (2003) in a joint effort of the Univ. 
of Chicago group headed by Jon Rosner, along 
with one of the current authors (P. Gorham) of 
this chapter. This study investigated a 10 19 eV 
vertical air shower, including explicit geomag- 
netic effects, with general interest in elucidating 
issues for detection by a possible radio augmen- 
tation to the Auger Observatory for ultra-high 
energy cosmic rays. 

The Suprun et al. simulation did not make any 
simplifying assumptions regarding far-field con- 
ditions. Instead, the electrodynamics simulation 
began with the general formula for a radiating 
particle (Jackson 1999, Zas et al. 1992) in arbi- 
trary motion: 



n x (n — n(3) x (3 



|l-n/3-n| 3 Z 



(4) 



E(x,i a ) 



n/3 



7 2 |l-n/3-n| 3 Z 2 



which is correct regardless of the distance to the 
antenna. In this formula f3 is the velocity vector 
in units of c, (3 = d/3/dt is the acceleration vector, 
divided by c, n is a unit vector from the radiating 
particle to the antenna, and I is the distance to 
the particle, /i ~ 1 denotes the relative magnetic 
permeability of air, n the index of refraction. The 
square brackets with subscript "ret" indicate that 
the quantities in the brackets are evaluated at the 
retarded time, not at the time t a when the signal 
arrives at the antenna. 

The first term decreases with distance as l/l 2 
and represents a boosted Coulomb field. It does 
not produce any radiation. The magnitudes of 
the two terms in Eq. (0} are related as l/(j 2 l) 
and \0\/c. The characteristic acceleration of a 
30 MeV electron (7 w 60) of an air shower in 
the Earth's magnetic field (B « 0.5 Gauss) is 
|a| = ecB/{^m) ss 4.4 • 10 13 m/s 2 . Even when 
an electron is as close to the antenna as 100 m, 
the first term is two orders of magnitude smaller 
than the second and can be neglected. The sec- 
ond term falls as l/l and is associated with a ra- 
diation field. It describes the electric field of a 
single radiating particle for most geometries rel- 
evant to extensive air showers. It can be shown 
(Wheeler & Feynmann 1949) to be proportional 
to the apparent angular acceleration of the charge 
up to some non-radiative terms that are propor- 
tional to l/l 2 . This relation is referred to in the 
literature as "Feynman's formula." 

Suprun et al. did not, however, yet perform 
a full cascade calculation, but rather used a 
parametrisation of the shower density to gener- 
ate a shower profile, then used Monte Carlo tech- 
niques to sample the particle distribution obey- 
ing this parametrisation. In one of the longest- 
standing empirical models for air shower devel- 
opment, called the Nishima, Kamata, Greisen 
(NKG) model, the lateral particle density p e 
is parametrized by the age parameter s of the 
shower (s = 1 for the shower maximum) and the 
Molicre radius r m (Bourdeau et al. 1980, Greisen 
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1956, Kamata & Nishimura 1958): 

/ \ s— 2 / \ s— 4.5 

Pe=K N ( ^— 1 + , (5) 

\s m T m j \ s 7n r m j 

where 

„ TV r(4.5-a) 

w 27r«r( s )r(4.5-2 S )' lj 

r is the gamma function, r the distance from 
the shower axis, TV the total number of charged 
particles, and s m = 0.78 — 0.21s. The Moliere 
radius for air is approximately given by r m = 
74 (po/ p) m, with po an d p being the air densities 
at sea level and the altitude under consideration, 
respectively. 

As a shower travels toward the Earth and en- 
ters denser layers of the atmosphere, the age pa- 
rameter increases while the Moliere radius drops. 
Both processes affect the spread of the lateral dis- 
tribution. The influence of the age parameter ap- 
pears to be more significant. As it grows, the 
average distance of the shower particles from its 
axis increases. This effect overcomes the influ- 
ence of a smaller Moliere radius which tends to 
make the lateral distribution more concentrated 
toward the axis. For a fixed age parameter s, 
however, the Moliere radius is the only quantity 
that determines the spread of the lateral distri- 
bution. At shower maximum (s = 1) the average 
distance from the axis can be calculated to be 
(2/3)s m r m = 0.38 r m . 

Fig. [§] shows the results of the Suprun et al. 
simulation for the shape of the intrinsic radio 
pulse, in terms of field strength vs. time at the 
receiving antenna location, though without any 
of the filtering effects of any antenna imposed on 
it yet. Fig. EH gives the Fourier transform E v 
of this pulse. The nonzero thickness of the air- 
shower pancake translates into a loss of coherence 
at frequencies corresponding to wavelengths com- 
parable to the shower thickness, thereby limiting 
the main part of the radiation spectrum to the 
frequencies below 100 MHz. 

These simulations, though using a greatly 
thinned set of input particles (10 4 compared to 
10 10 in actuality) do show characteristics similar 
to what was observed historically. In addition, 
the simulations also begin to reveal some of the 
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Figure 9. The EW component Eew of electro- 
magnetic pulse of 0.33 • 10 10 radiating electron- 
positron pairs distributed over the thickness of 
the shower pancake at 1800 m above sea level. 
The axis of the pancake is located 200 m South 
of the antenna. The time axis was chosen in such 
a way that the pulse produced by a pair located 
in the axis at the bottom of the pancake starts at 
time 0. 



geomagnetic complexity of the emission pattern, 
suggesting reasons for some of the surprising vari- 
ations observed in the measurements by antenna 
arrays. 

Consider the frame centered at the antenna, 
with axis Ox going to the magnetic West, Oy to 
the South and Oz directly up. The initial velocity 
of all charged particles is assumed to be vertical: 
(3 = (0,0,-1), while the initial acceleration f3 is 
parallel to Ox, or, in other words, to the (1, 0, 0) 
vector. 

Electrons bend toward the magnetic West and 
positrons toward the East. The electric fields 
from both particles of an electron-positron pair 
are coherent; the opposite signs of their acceler- 
ations are canceled by the opposite signs of the 
electric charges. 
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Figure 10. The Fourier transform of the electro- 
magnetic pulse shown in Fig. El The spectrum is 
very flat below 2 MHz. The limited statistics of 
the model results in some jitters at 200—500 MHz. 
The spectrum above 500 MHz is not shown be- 
cause the statistics is not sufficient to make re- 
liable predictions of the Fourier components at 
these high frequencies. 



Let ip be the angle between Ox and the direc- 
tion to the shower core, R the distance to the 
core, and h the altitude of the radiating parti- 
cle above the antenna. The denominator of the 
second term of Eq. (@J is independent of ip. The 
numerator determines that, to leading (second) 
order in R/h, the initial electric field vector E re- 
ceived at the antenna lies in the horizontal plane 
and is parallel to (cos 2ip, sin 2ip, 0) (Green et al. 
2003): 

E || (cos 2<0, sin 2-0,0). (7) 

The magnitude of the numerator is independent 
of the angle ip up to terms of order R A /h 4 . This 
result shows that although particles are acceler- 
ated by the Earth's magnetic field in the EW di- 
rection regardless of angle ip, the radiation re- 
ceived at the antenna does not show preference 
for the EW polarisation. Instead, it is directly 
related to the angle ip. As the particle trajectory 



bends in the Earth's magnetic field and the veloc- 
ity deflects from the vertical direction, the rela- 
tion Q between the direction of the electric field 
vector and angle ip does not hold. Nonetheless, it 
will be useful for understanding the angular de- 
pendence of the electric field. 

Suprun et al. computed electromagnetic pulses 
for the pancakes with axes located at the same 
distance R — 200 m from the antenna but at 
various angles ip from the Ox direction. Fig. 1111 
shows the radio signal strengths that would be 
received by EW and NS-oriented antennas. Note 
that Eq. JJJ) predicts that components of the radi- 
ation coming from the start of the particle trajec- 
tory vanish at some angles ip: Eew = at ip — 
±71-/4, ±3tt/4, while E NS = at ip = 0, ±tt/2, tt. 
This fact explains why E u ew is relatively small at 
ip — ±7r/4, ±37r/4 and E u ms is small at ip = 0, tt 
(Fig. ITT)). Another mechanism is responsible for 
E v ns being virtually at ip — ±7r/2. At these 
angles the trajectories of two charged particles of 
an electron-positron pair are symmetric with re- 
spect to the yOz plane. The NS component of 
radiation emitted by this pair vanishes not only 
at the start but throughout its flight. 

Modeling by the LOPES collaboration. 

Another simulation effort is under way in the 
Max-Planck- Institut fr Radioastronomie at Bonn, 
led by T. Huege and an author of this chapter 
(H. Falcke). This group is part of a collaboration 
developing the LOFAR Prototype Experimental 
Station (LOPES), an engineering model of one 
station of the Low Frequency Array (LOFAR). 
LOPES is operating jointly with the KASCADE 
Grande air shower array in Karlsruhe (Schieler et 
al. 2003); LOFAR is a funded effort to develop a 
very large area ground array for radio astronomy 
in the HF to VHF regime, sharing many com- 
mon interests with the SKA. The LOPES group 
has recently published a detailed analysis of the 
geosynchrotron model for the case of a 10 17 eV air 
shower (Huege & Falcke 2002, 2003) in prepara- 
tion for a major effort at an electro-dynamical air 
shower Monte Carlo code (Huege & Falcke 2004). 

The LOPES group has taken special care of 
taking into account the longitudinal development 
of the air shower by performing an integration 
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Figure 1 1 . The East- West and North-South com- 
ponents of the field strength \E v ew\ and \E u ns\ 
(circles and triangles, respectively) at 55 MHz 
as functions of angle ip between the magnetic 
West and direction to the shower core. The dis- 
tance between the origin and a circle or a tri- 
angle represents the field strength in the units of 
/zV/m/MHz. The angular spacing between circles 
or triangles is n/8. At V = ±7r/2 \E v ^s\ do not 
exceed 0.1 //V/m/MHz and two triangles overlap. 
All points were calculated for the vertical shower 
at a 200 m distance from the antenna. 



over the shower as a whole, and they have con- 
sidered the variation of the field strength as a 
function of radial distance from the shower core 
as well. They use a shower parametrisation based 
on the NKG model with a shower disk that flares 
out from the center, in a manner similar to the 
Chicago/Hawaii study, and thus, apart from the 
energy difference, the results do bear some com- 
parison. The LOPES study also did an integral 
over a power-law distribution of electron energies, 
appropriate to an air shower. However, they did 
not do any near-field corrections to their results, 



but this is not a major drawback for a lower en- 
ergy shower since these showers do reach their 
maxima at altitudes of typically several km away 
from an observer on the ground. 

Fig. 1121 shows the spectrum emitted by the air 
shower maximum for a shower disk profile with 
realistic flaring according to the parametrisations 
of Agnetta et al. (1997) and Linsley (1986). As 
expected, the spectrum emitted by the Linsley 
flaring disk extends to higher frequencies than the 
one generated by the Agnetta flaring disk because 
of the lower thickness in the shower centre where 
most of the particles reside. 
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Figure 12. 



E(R, 2irv) 



-spectrum at the centre 



of the area illuminated by the maximum of a 
10 17 eV air shower with flaring T-pdf, i?o = 4 km 
and a broken power-law energy distribution from 
7 = 5-1000. Solid: flaring (Agnetta et al. 1997) 
lateral distribution, short-dashed: flaring (Lins- 
ley 1986) lateral distribution 



The modeled radial dependence at different fre- 
quencies is shown in Figure 1131 Here the three 
families of curves represent different frequencies, 
and the different slopes between the two curves 
at a given frequency are for the cases of an ob- 
server with a given distance from the shower cen- 
ter in the directions perpendicular and parallel 
to the geomagnetic field. This result thus indi- 
cates again the importance of the geomagnetic 
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effects in the azimuthal distribution of radiation 
for a given magnetic field direction. Early results 
from the upcoming detailed Monte Carlo simula- 
tions of the LOPES collaboration, however, show 
that asymmetries in the emission pattern due to 
the geomagnetic field seem to be washed out to a 
high degree once realistic distributions of particle 
track lengths are taken into account. 
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Figure 14. Reconstructed pulses emitted by the 
maximum of a 10 17 eV shower with flaring (Ag- 
netta et al. 1997) T-pdf, broken power-law energy 
distribution from 7 = 5-1000 and Rq — 4 km, 
using an idealized rectangle filter spanning 40- 
160 MHz. Solid: centre of illuminated area, short- 
dashed: 100 m to north from centre, dash-dotted: 
250 m to north from centre 



Figure 13. Radial dependence of E(R, 2nv) for 

the maximum of a 10 17 eV air shower with flaring 
(Agnetta et al. 1997) T-pdf, R = 4 km and a 
broken power-law energy distribution from 7 = 
5-1000. Solid: v = 50 Mhz, short-dashed: v = 
75 Mhz, long-dashed: v — 100 Mhz, upper curves 
for distance from shower center to the east-west, 
lower curves for distance to north-south. 



Fig. El shows a reconstructed pulse generated 
by the flaring Agnetta disk as it would be mea- 
sured by a receiver with a given bandwidth. The 
pulse amplitude drops noticeably when the ob- 
server moves from the centre of the illuminated 
area on the ground to a distance of 100 m, and is 
already quite diminished at a distance of 250 m. 

The LOPES study addresses the important 
problem of integrating over the shower evolu- 
tion as a whole in a simplified fashion by ap- 
proximating the shower evolution with a num- 
ber of discrete steps. The characteristic scale 



for these steps is given by the "radiation length" 
of the electromagnetic cascades in air, Xq = 
36.7 g cm -2 , corresponding to « 450 m at a 
height of 4 km. One can therefore discretise the 
shower evolution into "slices" of thickness Xq, as- 
suming these contain independent generations of 
particles and therefore radiate independently. Su- 
perposition of the individual slice emissions, cor- 
rectly taking into account the phases arising from 
arrival time differences, then leads to the total 
emission of the shower. 

For a vertical 10 17 eV air shower at a height of 
i?o = 4 km they add the emission from eight slices 
above and eight slices below the shower maximum 
to the emission from the maximum itself. The 
closest slice then lies at i?o = 950 m from the 
observer, a distance they did not want to fall be- 
low because of approximations contained in their 
calculations that are only valid in the far-field. 

Although this treatment is clearly oversimpli- 
fied, the results depicted in Fig. El indicate that 
the integration over the shower as a whole signif- 
icantly enhances the emission strength and thus 
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cannot be neglected. In particular, this implies 
that the emission is actually not dominated by 
a narrow region around the shower maximum, 
but that the entire shower evolution contributes. 
A realistic treatment of the integration over the 
shower as a whole is carried out as part of the up- 
coming Monte Carlo simulations of the LOPES 
collaboration (Huege & Falcke 2004). 

Data from the LOPES Experiment will become 
available soon, but first results indeed confirm 
the association of the air shower with a sharp ra- 
dio pulse, having the expected properties (e.g., 
Horneffcr ct al. 2004). This puts the radio detec- 
tion method on rather firm ground. 
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Figure 15. E(R, 2nv) -spectrum of a full (lon- 
gitudinally integrated) 10 17 eV air shower with 
flaring (Agnetta et al. 1997) T-pdf, R — 4 km 
and a broken power-law energy distribution from 
7 = 5-1000. Solid: centre of illuminated area, 
short-dashed: 100 m to north from centre, long- 
dashed: 250 m to north from centre, black points: 
re-scaled (Spencer 1969) data as presented by (Al- 
lan 1971), grey points: re-scaled (Prah 1971) data 



It is interesting that in spite of the differences in 
the approach from the LOPES studies and those 
of the Chicago/Hawaii group, the results for the 
radio spectrum for a distance of 200/250 m from 
the shower core show a very similar frequency de- 



pendence, with the field strength falling about a 
factor of 300 as one goes from 10 to 100 MHz. 
The absolute value of the field strength is about 
a factor of 30 or so different, which is inconsistent 
with a strict linear scaling of field strength with 
energy as one might expect; however, agreement 
to within a factor of 2-3 is actually quite good 
considering the fact that these are completely in- 
dependent efforts. 

2.3. Askaryan effect and its confirmation 

As noted early in this discussion, the Askaryan 
effect was the original motivation for much of the 
effort to measure radio emission from air show- 
ers, but the coherent geo-synchrotron emission 
detailed above was found to be the dominant 
contribution for air showers, and the coherent 
Cherenkov emission from the charge excess, while 
not discounted, was largely forgotten because of 
its small contribution. However, for showers in 
solid materials such as ice or the lunar regolith 
which are relatively radio-transparent, the shower 
lengths are short enough (~ 10 m) that the mag- 
netic effects leading to synchrotron emission may 
be neglected, and the coherent Cherenkov emis- 
sion becomes the more important secondary radi- 
ation. Here the quadratic rise of radio power with 
frequency leads to the conclusion that, at energies 
above 10 18 eV, the coherent Cherenkov emission 
will dominate all secondary radiation, including 
optical emission, by a wide margin. 

Although it is not presently possible to produce 
EeV cascades in terrestrial accelerators, electro- 
magnetic showers with composite total energies 
in this range can be easily synthesized by super- 
posing gamma-rays of energies above the pair- 
production threshold. If the gamma-ray bunch 
is small compared to the wavelength of the radio 
emission (true for most pulsed linacs), the result- 
ing showers will differ from natural EeV showers 
only logarithmically, due to the details of the ini- 
tial interaction. However, since the bulk of the 
radio emission arises from the region of maxi- 
mum shower development, the differences in radio 
Cherenkov emission are modest and easily quan- 
tified. 

In mid-2000, Askaryan's hypothesis was in fact 
confirmed at the Stanford Linear Accelerator 
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Center (SLAC) in an experiment using a silica- 
sand target and pulsed gamma-ray bunches with 
composite energies in the EeV range (Saltzberg 
et al. 2001). In the 2002 follow-on experiment 
(Gorham et al. 2004), the sand was replaced by 
synthetic rock salt, which has a higher dielectric 
constant and lower loss tangent than silica sand, 
and further studies were made of the polarisation 
behavior of the emission. 

Fig. shows a typical pulse profile (in- 

set) and a set of measured peak field strengths 
for pulses taken at different points along the 
shower in the 2000 experiment. The plot- 
ted curve shows the expected profile of the to- 
tal number of particles in the shower, based 
on the Kamata-Nishimura-Greisen approxima- 
tion (Saltzberg 2001). Here the field strengths 
have been scaled in the plot to provide an ap- 
proximate overlay to the relative shower profile. 
Clearly the pulse strengths are highly correlated 
to the particle number profile. Since the ex- 
cess charge is also expected to closely follow the 
shower profile, this result confirms Askaryan's hy- 
pothesis. 

Pulse polarisation was measured with an S- 
band (2 GHz) horn directed at a shower posi- 
tion 0.5 m past the shower maximum. Fig. I16B 
shows the pulse profile for both the 0° and 90° 
(cross-polarised) orientations of the horn. The 
lower two panes of this portion show the derived 
degree of linear polarisation and the angle of the 
plane of polarisation, respectively. Because of the 
vector correlation of the pulse polarisation with 
the shower velocity vector and the Poynting flux 
vector, it is possible to use the angle of the polar- 
isation to track the shower axis. An example of 
this is shown in Fig.llGC , where the angle of the 
plane of polarisation is plotted at three locations 
with respect to the shower axis, showing the high 
correlation with the predicted angle. 

Fig.lltjU shows a typical sequence of pulse field 
strengths versus the total shower energy. The fit- 
ted linear rise of field strength with beam current 
is consistent with complete coherence of the radi- 
ation, implying the characteristic quadratic rise 
in the corresponding pulse power with shower en- 
ergy. Fig.llfii'' shows a similar result for the 2002 
experiment, but now covering a much wider range 



of energy, plotted as pulse power instead of field 
strength. The Askaryan process is found to be 
quadratic over four orders of magnitude in shower 
energy. 

Fig. I16E shows the spectral dependence of 
the radiation, which is consistent with the lin- 
ear rise with frequency that is also characteristic 
of Cherenkov radiation. Also shown is a curve 
based on a parametrisation of Monte Carlo re- 
sults (Zas et al. 1992). The uncertainties are es- 
timates of the combined systematic and statisti- 
cal uncertainties. Note that the figure compares 
absolute field strength measurements to the pre- 
dictions and the agreement is very good. 

In summary, there is clear experimental ev- 
idence that Askaryan's hypothesis is confirmed 
and that the predicted emission from high energy 
cascades is present in the expected amounts. This 
lends strong support to experiments designed to 
exploit this effect for high energy neutrino and 
cosmic ray detection. 

3. Prospects for the SKA 

The ultra-fast transient radio events described 
here will either be signal, or at some level, back- 
ground for the SKA. With recent pulsar studies 
extending to broader and broader bandwidths, 
and faster and faster pulse transients, under- 
standing of these events may become important in 
verifying the detection of pulsar transient events, 
certainly a mainstay of SKA scientific interest. 
Whether or not the SKA can be a competitive 
instrument for the detection of the various types 
of events described here will depend strongly on 
the final choice of design. However, without care- 
ful choices made at this early stage of the effort, it 
will be much more likely that the design is "pes- 
simized" rather than optimized for their detec- 
tion. 

3.1. Cosmic ray air shower detection 

From the description above, it is evident that 
air shower radio emission in the > 10 18 eV energy 
regime has three important characteristics which 
will impact the design of any radio array with 
intent to detect them: 

• The optimal frequency range is ~ 20 — 
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Figure 16. A: Shower RF field strength profile with typical pulse (inset). B: Polarisation measurements 
of a typical RF Cherenkov pulse at 2 GHz. C: Correlation of plane of polarisation with antenna offset 
from shower axis. D: Coherence of RF Cherenkov at 2 GHz, measured during 2000 SLAC experiment. E: 
absolute field strength and prediction from Cherenkov. F: Coherence of radiated power over the 0.2-1.2 
GHz band. 
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200 MHz; 

• The physical area over which one expects to 
detect the emission is limited in diameter to 
a few tens of km, and often much less. 

• The time-scale for the air shower radio emis- 
sion is an impulse of order 20 ns or less in 
duration. 

An SKA design extending down to 200 MHz 
is still adequate, though not optimized, for air 
shower detection, but the sparsity of the array 
and the bandwidth of the front-end receivers will 
have significant impact on the triggering and re- 
construction of the shower energy and direction, 
and of course the ultimate sensitivity. For the 
SKA to be competitive for giant air shower detec- 
tion, careful consideration of all of these factors 
will make the difference as to whether the SKA is 
irrelevant for this field, or a dominant competitor. 

Because of the wide variety of designs currently 
under consideration for the SKA, and their rapid 
evolution, it is impractical to assess each one for 
its capability in air shower radio detection. In- 
stead, we take the approach of estimating what 
would make the ideal detector for air shower ra- 
dio detection, and then we consider how this com- 
pares to current plans. 

As noted above, one would ideally like to work 
at frequencies that extend below 100 MHz, but 
with as broad a bandwidth as possible, able even 
to resolve the ~ 10 ns time-scale for air shower 
radio emission. This requires several hundred 
MHz of bandwidth, extending perhaps down to 
50 MHz or so. The immediate implication is that 
the fundamental detector element must be a very 
broad-band antenna, perhaps with of order 6:1 
bandwidth or more. A high-gain antenna is also 
problematic, since one cannot predict a priori the 
arrival direction of an air shower radio pulse. For 
that reason a phased-array concept with digital 
beam-forming is by far superior to other designs. 
Dual polarisation is also clearly desirable, since 
the radiation itself is highly polarised. 

These considerations lead one to consider scale- 
invariant designs for the primary antenna element 
such as dual-polarisation log-periodic dipolc an- 
tennas (LPDAs), but their typical beam- widths 



(60 — 90°) would require a cluster of at least 6 
antennas to get even minimal coverage of the 
entire sky. LPDAs have one other character- 
istic which is undesirable for air shower detec- 
tion: their inherent pulse dispersion reduces sen- 
sitivity to impulsive events unless a de-dispersing 
compensator either analog or digital) is imple- 
mented. However, if this can be overcome, they 
are lightweight, easy to construct and straight- 
forward for impedance matching and modeling. 
Variations on the LPDA design could also be 
scaled up as stand-alone units to satisfy this need, 
and are more compact with potentially better 
phase centers. A non-dispersive alternative could 
be an array of quad-ridged horns, which can rou- 
tinely achieve the 6:1 bandwidth required, but 
they would be larger and heavier than LPDAs. 
Obviously, even a simple active "inverted- V" an- 
tenna as used for LOFAR and LOPES are also 
very useful if optimized for the right frequencies. 

The use of such a broad-band system of course 
raises the question of how one can possibly deal 
with interference. This has been successfully 
demonstrated with the LOPES experiment (see 
Horneffer et al. 2003). Another excellent exam- 
ple of a solution to this problem is the FORTE 
satellite (Lchtincn et al. 2004), which launched 
in 1997 with a 30-300 MHz nadir-pointing dual- 
polarisation LPDA with a tunable 25 MHz receiv- 
ing band. FORTE was optimized for detection of 
electromagnetic impulsive events, and its mission 
was to provide an unclassified test-bed for nuclear 
treaty verification efforts while pursuing a science 
program of lightning and atmospherics detection. 

At an orbital altitude of 800 km, FORTE was 
constantly exposed to a barrage of anthropogenic 
EM interference. FORTE was able to retain trig- 
gering capability for impulses down to a level 
within about 5cr of the ambient thermal noise 
level by sub-dividing their large band into a se- 
ries of 1 MHz channels and triggering when a 
majority of the bands exceeded threshold indi- 
cating a broad-band pulse. The signal digitiza- 
tion was still done over the entire bandwidth, pre- 
serving the broad-band coherence of the impulse. 
But since the vast majority of anthropogenic in- 
terference is inherently narrow-band, the multi- 
band trigger technique was very effective, when 
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combined with a so called noise-riding threshold 
which effectively maintained the trigger rate for 
each sub-band to a constant level. This greatly 
reduced the ability of strong narrow-band carriers 
to cause rapid re-triggering of one of the channels 
which might skew the broad-band trigger rate. 
As a result, analysis of FORTE data has recently 
even provided the first published limits on neu- 
trino fluxes in energy regimes of <~ io 22 ~ 24 eV, 
based on the lack of observed radio impulses 
emanating from within the Greenland ice sheet 
(Lehtinen et al. 2004). 

Applied to a potentially much broader-band 
system as proposed above for SKA air shower 
radio detection, the multi-band triggering would 
in principle be applied to each cluster locally. 
If a trigger occurred, it would cause a global 
broadcast out to stations within a several km ra- 
dius of the triggered cluster, interrogating these 
other stations to see if they also triggered., When 
enough stations trigger to justify it, a global trig- 
ger would be initiated and all of the stations 
within the affected distance (including appropri- 
ate margin to establish the boundaries of the af- 
fected area) would save their buffered data. 

The design implications for such a system 
clearly favor the phase aperture array concept for 
the SKA, with a low frequency cutoff at the lower 
end of the VHF band. In many ways the con- 
cept of a Low Frequency Array (LOFAR) is per- 
haps best matched to air shower radio detection, 
and insight can be gained toward adaptation of 
the SKA possibilities by considering the adapta- 
tion required for a LOFAR-type array. Studies 
for such applications with LOFAR have been re- 
cently published (Falcke & Gorham 2003, Huege 
& Falcke 2002, 2003), and the results are quite 
promising. 

3.2. Neutrino detection 

In contrast to the problem of air shower radio 
detection with the SKA, which is driven by the 
fact that there is only one clear mechanism for de- 
tection, neutrino detection with the SKA may be 
pursued on several fronts. The scientific motiva- 
tions for both neutrino and air shower detection 
from EeV to ZeV energies are closely related, and 
neutrino detection at these energies will provide 



highly complementary information to our current 
incomplete knowledge of the sources and propa- 
gation of the highest energy cosmic rays. 

To date, no cosmic high energy (> 1 GeV) 
neutrinos have been detected from any source. 
The AMANDA detector at Amundsen station, 
Antarctica, has detected cosmic-ray secondary 
neutrinos up to ~ 100 TcV energies, but these 
arise from interactions of garden-variety ~ PcV 
cosmic rays in the Earth's atmosphere. 

For this reason, the discussion of neutrino de- 
tection must be more broad in scope, since we do 
not yet know which detection channels might lead 
to methods with sufficient sensitivity to see fluxes 
of neutrinos over the entire range of 10 10-23 eV 
where they are expected but so far unobserved. 
This section is therefore more speculative with 
regard to possible techniques, but appropriate to 
the high level of scientific interest in neutrino de- 
tection. 

3.2.1. Neutrino interactions in the Earth 

At energies of about 1 PeV, the earth becomes 
opaque to neutrinos at the nadir. For higher en- 
ergies, the angular region of opacity grows from 
around the nadir till at EeV energies, neutrinos 
can only arrive from within a few degrees below 
the horizon. The interaction length at these en- 
ergies is of order 1000 km in water, so such neu- 
trinos have a significant probability of interacting 
along a ~ 100 km chord. If the interaction takes 
place within ~ 10 m of the surface of rock or dry 
sand or soil, the resulting cascade will produce 
coherent Cherenkov radiation up to microwave 
frequencies. Thus, for example, since arrays are 
often sited with mountains or ridges near the hori- 
zon, the entire near-surface volume of the moun- 
tain range becomes a neutrino target, and events 
can originate anywhere along its surface. The flux 
density expected for such events (cf. Saltzberg et 
al. 2001) is 

S v * 12 MJy (^-) 2 x 

(lO 18 ev) (2OO MHz) ^ 
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where E c is the cascade energy and R the dis- 
tance to the cascade. 2 The Cherenkov process 
weights these events strongly toward the higher 
frequencies, though events that originate deeper 
in the ground will have their spectrum flattened 
by the typical v~ Y behavior of the loss tangent of 
the material. 

A similar process leads to coherent transition 
radiation (TR; cf. Takahashi et al. 1994) from 
the charge excess of the shower, if the cascade 
breaks through the local surface. TR has spec- 
tral properties that make it more favorable for 
an array at lower frequencies: it produces equal 
power per unit bandwidth across the coherence 
region. The resulting flux density for a neutrino 
cascade breaking the surface near the array ar- 
ray, observed at an angle of within ~ 10° from 
the cascade axis, is (cf. Gorham et al. 2000): 

s v ,tr{o< io°) ~ 2MJy (riL) x 

The implication here is that, if an array can retain 
some response from the antennas to near-horizon 
fluxes, the payoff may be a significant sensitivity 
to neutrino events in an energy regime of great in- 
terest around 1 EeV, or even significantly below 
this energy depending on the method of trigger- 
ing. 

3.2.2. Neutrino interactions in the atmo- 
sphere 

Neutrinos can themselves also produce air 
showers. The primary difference between these 
and cosmic-ray-induced air showers is that their 
origin, or first-interaction point, can be anywhere 
in the air column, with an equal probability of 
interaction at any column depth. Neutrino air 

2 Note that in this case the neutrino energy is not necessar- 
ily equal to the cascade energy E c , because for the typical 
deep-inelastic scattering interactions that occur for EeV 
neutrinos, only about 20% of the energy is put into the 
cascade, while the balance is carried off by a lepton. For 
electron neutrinos, the electron will rapidly interact and 
add its energy to the shower, but for muon or tau neutri- 
nos, this lepton will generally escape undetected (although 
the tau lepton will itself decay within a few tens of km at 
1 EeV). 



showers can even be locally up-going at mod- 
est angles, subject to the earth-shadowing effects 
mentioned above. 

Detection of such events is identical to detec- 
tion of cosmic-ray-induced air showers, except for 
the fact that sensitivity to events from near the 
horizon is desirable, since these will be most eas- 
ily distinguished from cosmic-ray-induced events. 
Beyond a zenith angle of ~ 70° cosmic-ray ra- 
dio events will be more rare, and those that are 
detected in radio will be distant. The column 
depth of the atmosphere rises by a factor of 30 
from zenith to horizon; thus cosmic ray induced 
air showers have their maxima many kilometers 
away at high zenith angles. Neutrino showers in 
contrast may appear close by, even at large zenith 
angles. 

Of particular interest is the possibility of ob- 
serving "double-bang" (Learned & Pakvasa 1995) 
tau neutrino events. In these events, a v T inter- 
acts first, producing a near-horizontal air shower 
from a deep-inelastic hadronic scattering interac- 
tion. The tau lepton escapes with of order 80% 
of the neutrino energy, and then propagates an 
average distance of 50£? T /(10 18 eV) km before 
decaying and producing (in most cases) another 
shower of comparable energy to the first. Detec- 
tion of both cascades within the boundaries of a 
surface radio array would provide a unique sig- 
nature of such events. And in light of the recent 
neutrino results indicating — > v T oscillations, 
it is likely that neutrinos from astrophysically dis- 
tant sources would be maximally mixed, leading 
to a significant rate of v T events. 

3.2.3. Neutrino interactions in the lunar 
regolith 

There is an analogous process to the earth- 
surface layer cascades mentioned above which can 
take place in the lunar surface material (the re- 
golith). In this case the cascade takes place as 
the neutrino nears its exit point on the moon af- 
ter having traversed a chord through the lunar 
limb. This process, first suggested by Dagke- 
samansky & Zheleznykh (1989) is the basis of 
several searches for diffuse neutrino fluxes at en- 
ergies of - 10 20 cV (Hankins et al. 1996; Gorham 
et al. 1999, 2001, Gorham et al. 2003) using 
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large radio telescopes at microwave frequencies. 
Based on the simulations for these experiments 
(Alvarez Muniz & Zas 1997a, 1997b; Zas, Halzen 
& Stanev 1992) and confirmation through several 
accelerator measurements (Gorham et al. 2000 ; 
Saltzberg et al. 2001), the expected flux density 
from such an event at about 1 attenuation-length 
depth in the regolith can be roughly estimated as 




Note here that the flux density is far lower than 
for air shower events, but the two should not be 
compared, since the lunar regolith events are co- 
herent over ~ degree angular scales, correspond- 
ing to several thousand km at the Earth's surface. 
They also originate from a small, known angular 
region of the sky (the surface of the moon) . Thus 
their detectability depends on the sensitivity of 
the synthesized beam, and on the ability of the 
system to trigger on band-limited pulses. 

Transition radiation events may also be de- 
tectable in a similar manner, as noted above. For 
TR from events that break the lunar surface, the 
resulting pulse differs from a Cherenkov pulse be- 
cause it is flat-spectrum. Because TR is strongly 
forward beamed compared to the Cherenkov radi- 
ation from the moon, we estimate that the max- 
imum flux density for this case, at an angle of 
<~ 1.5° from the cascade axis, is about a factor 
of 20 higher than at ~ 10°. At earth the implied 
flux density for LOFAR is: 

S max , TR (6^1.5°) ~ 40 Jy (j^-^j -(H) 

Although this channel does not provide a higher 
flux density than the Cherenkov process, it is a 
flat spectrum process that may in some cases pro- 
vide more integrated flux across a given band. 

These pulses are essentially completely band- 
limited prior to their entry into the ionosphere, 
with intrinsic width of order 0.2 ns. Dispersion 
delay in the ionosphere will of course significantly 
impact the shape of any pulse of lunar origin. 
This will limit the coherence bandwidth for a 
VHF system. The dominant quadratic part of 



the dispersion gives an overall delay 

r lon = 1.34 xHT 7 ^ (12) 

where Ti on is the delay in seconds at frequency 
v (in Hz) for ionospheric column density N e in 
electrons per m 2 . For typical nighttime values of 
N e ~ 10 17 m~ 2 the zenith delay at 200 MHz is 
330 ns, and the differential dispersion is of order 
3 ns per MHz, increasing at lower frequencies as 
. For bandwidths up to even several tens of 
MHz for zenith observations, and perhaps a few 
MHz at low elevations, the pulses should remain 
band-limited. However, coherent de-dispersion 
will be necessary to accurately reconstruct the 
broad-band pulse structure. 

Although the problem of coherent de- 
dispersion is a difficult one, a system operating in 
the 0.2-1 GHz range may have an edge in sensitiv- 
ity over systems operating at higher frequencies, 
under conditions where the intrinsic neutrino 
spectra are very hard. This is due to the fact 
that the loss tangent of the lunar surface mate- 
rial is relatively constant with frequency (Olhoeft 
and Strangway 1976), and thus the attenuation 
length increases inversely with frequency. This 
means that a lower frequency array may probe 
a much larger effective volume of mass than the 
higher frequencies can. At 200 MHz, the RF 
attenuation length should be of order 50 m or 
more, compared to 5-7 m at 2 GHz. When this 
larger effective volume is coupled with the larger 
acceptance solid angle afforded by the broader 
RF beam of the low-frequency Cherenkov emis- 
sion, the net improvement in neutrino aperture 
could well compensate the loss of sensitivity at 
lower energies by a large margin. 

It is also worth noting here that these lunar re- 
golith observations are distinct from other meth- 
ods in high energy particle detection, in that they 
do require the array to track an astronomical tar- 
get, and can and will make use of the synthetic 
beam of the entire array. This is because, al- 
though the sub-array elements should be used for 
the detection since they will have a beam that 
covers the entire moon, the Cherenkov beam pat- 
tern from an event of lunar origin covers an area 
of several thousand km wide at earth, and is thus 
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broad enough to trigger the entire array. Post- 
analysis of such events can then localize them 
to a few km at 100 MHz and about 200 m at 
1.4 GHz on the surface of the moon, providing 
opportunities for more detailed reconstruction of 
the event geometry. If, as expected, a high res- 
olution (5 m) 3D cartographic map of the moon 
will be produced by the Terrain Mapping Cam- 
era of India's Chandrayaan-1 lunar mission due 
for launch in 2008 [www.isro.org/chandrayaan- 
1/announcement.htm] becomes available, this 
may be used to determine the local gradient and 
roughness of the surface near the position of cas- 
cade exit. We anticipate that this information, 
together with the polarisation angle and its fre- 
quency dependence, both of which can be mea- 
sured by the SKA, may enable neutrino direction 
reconstruction to an "event arc" on the sky of 
thickness a few degrees to be routinely achieved 
on an event-by-event basis. 

It is in the neutrino energy range most sensitive 
to UHECR origin at 10 9 -10 13 GeV that the SKA 
may have the greatest impact by using the lunar 
Cherenkov technique. So far no UHE neutrinos 
have been detected and the current observational 
limits are shown in Fig. 1171 

In the 10 n -10 13 GeV range the Goldstone Lu- 
nar Ultra-high Energy (GLUE) neutrino exper- 
iment (Gorham et al. 2003) has the best limit. 
Other planned experiments in this energy range 
such as SALSA (Gorham et al. 2002) and ANITA 
(Barwick et al. 2003) will lower these limits and 
hopefully detect neutrinos. The lunar Cherenkov 
technique used in the GLUE experiment was pi- 
oneered by Hankins, Ekers & O'Sullivan (1996) 
using the Parkes 64m radio telescope. The GLUE 
experiment used two dishes of the Goldstone 
Deep Space Tracking Network for 120 hours to 
look for Cherenkov radio emission from neutrino- 
induced cascades in lunar regolith. 

By scaling relationships given by Gorham et al. 
(2000) and Alvarez-Muniz & Zas (2001) describ- 
ing the electric field strength at the radio tele- 
scopes expected for a given cascade energy de- 
posited in the regolith (see also Beresnyak 2003), 
and comparing the proposed technical specifica- 
tions of the SKA (assuming 1 GHz frequency will 
be used, its higher bandwidth and larger tele- 
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Figure 17. Black solid curves with points show 
existing limits on diffuse neutrino intensity (for 
references see Gorham et al. 2003, from which this 
plot is adapted). Thick, dark red dotted curves 
show the expected sensitivity of Auger Observa- 
tory to and v e (top), and v T (with no deep- 
inelastic-scattering losses, bottom) in 5 years of 
observations. The magenta dash-dot curve la- 
belled SKA gives an estimates of the expected 
sensitivity of the SKA as a neutrino observato- 
ries for 120 h of lunar observations. ANITA is an 
Antarctic Long duration balloon mission due to 
fly in 2006; this estimate is for the full planned 3 
flights by 2009. Other curves bracket AGN pre- 
dictions, two TD models, and an estimate of the 
maximal GZK neutrino flux. 



scope field of view and collecting area) with those 
of the telescopes used in the GLUE experiment, 
one expects that for comparable lunar observing 
time with the SKA the threshold will be reduced 
to 2xl0 10 GeV and the sensitivity will be im- 
proved by a factor of about 2000 (dot-dash curve 
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"SKA" in Fig.ll7|l. making it potentially the most 
sensitive UHE neutrino observatory in the future 
for covering a large part of the important energy 
range 10 9 -10 14 GeV. 

4. Summary and outlook 

In the next several years giant air shower de- 
tectors will investigate the spectrum, composition 
and anisotropy of the UHECR, i.e. those with 
energies above 10 10 GeV, in an attempt to de- 
termine their origin. Because of its large collect- 
ing area, the use of the SKA to directly observe 
coherent geosynchrotron radio emission from cos- 
mic ray air showers has the potential to make a 
significant impact in this field. However, cosmic 
rays are deflected by magnetic fields and do not 
point back directly to their sources, and above 
~ 10 11 GeV UHECR suffer severe energy losses 
on interacting with CMBR photons, limiting their 
range to tens of Mpc from their sources. Hence 
studies of UHECR alone will probably be insuffi- 
cient to tie down their sources and whether they 
are accelerated or result from the decay of massive 
relic particles or emission by topological defects. 

UHE neutrinos are the key to determining the 
origin of these UHECR. This subject is of great 
importance to our understanding of the Universe 
as it impacts on our knowledge of dark matter, 
gravity, and high energy particle interactions. Di- 
rect radio observation by the SKA of air showers 
due to high energy neutrinos may contribute sig- 
nificantly to high energy neutrino astrophysics, 
particularly below 10 10 GeV. However, the enor- 
mous neutrino collecting area of the Moon, to- 
gether with the large aperture and excellent an- 
gular resolution of the SKA make UHE neutrino 
astrophysics using the lunar Cherenkov technique 
potentially the best approach for tying down the 
origin of the very highest energy particles in na- 
ture. 

The signal coincidence requirement between 
antennas and the nanosecond duration signal ex- 
perimental procedures are significantly different 
from those in normal radio astronomy, and must 
be taken into account together with the most ap- 
propriate signal processing technique for multi- 
ple antennas in the design of the SKA if it is to 



be used for lunar UHE neutrino observations and 
take a leading role in neutrino astronomy at the 
highest energies. 
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References 



[1] 
[2] 
[3] 



[G] 
[7] 



[9 

[10 
[11 

[12 
[13 

[14 
[15 



Abbasi, R. U., et al. 2004, Physical Review 
Letters, 92, 151101 

Agnetta, G. et al. 1997, Astropart. Physics, 
6, 301 

Allan, H.R., Neat, K.P., Jones, J.K. 1967, 
Nature, 215, 267 

Allan, H.R., Clay, R.W., Jones, J.K. 1970, 
Nature, 227, 1116 

Allan, H.R. 1971, Prog, in Elem. Part, and 
Cos. Ray Phys., ed. J. G. Wilson and S. A. 
Wouthuysen, (N. Holland Publ. Co.), Vol. 
10, 171 

Alvarez-Muniz, J., & Zas, E. 1997a, Phys. 



Lett. B, 411, 218; also astro-ph/9806098 



Alvarez-Muniz, J., & Zas, E. 1997b, Proc. 
XXVth ICRC, ed. M.S. Potgeiter et al. vol. 
7, 309 

Alvarez-Muniz, J., Zas, E. 2001, AIP Conf. 
Proc. 579, pp 128-138 

Alvarez-Muniz, J., Marques, E., Vazquez, 
R.A., Zas, E. 2003, Phys. Rev., D68, 043001 

Arons, J. 2003, ApJ, 589, 871 

Askar'yan, G.A. 1962, Sov. Phys. JETP, 14 
441; 1965, ibid 48, 988 

Barker,H.R., Neat, K.P., Jones, J.K. 1967, 
Nature 215, 267 

Barwick, S.W. et al. 2003, Proc. SPIE, 
4858, 265 



Beresnyak, A.R. 2003, astro-ph/0310295 
Bird, D.J. et al. 1995, ApJ, 441, 144 



25 



[16] Bourdeau, M.F., J.N.Capdevielle, J.N. & 
Procureur,J. 1980, J. Phys. G, 6, 901 

[17] Colgate, S.A., 1967, J. Geophys. Res., 72, 
4869 

[18] Dagkesamansky, R.D. & Zheleznykh, I.M. 
1989, Sov. Phys. JETP Lett., 50, 233 

[19] Dova, M.T., Fanchiotti, H., Garcia-Canal, 
C. et al. 1999, Proc. 26th ICRC, Salt Lake 
City Utah 1999, Ed. D. Kieda, M. Sala- 
mon, B Dingus, Vol. 1, p. 514 (HE.2.5.34) , 
http: / /krusty.physics.utah.edu/^icrcl999 



[20] Engel, R. Seckel, D. Stanev, T. 2001, Phys. 
Rev. D64, 093010 

[21] Falcke, H. & Gorham, P. 2003, Astropart. 
Phys., 19, 477 

[22] Gorham, P.W., Liewer, K.M., and Naudet, 
C.J. 1999, Proc. 26th ICRC, Salt Lake City, 
Utah, HE.6.3.15 

[23] Gorham, P.W, Saltzberg, D.P., Schoessow, 
P. et al. 2000, Phys. Rev., E 62, 8590 

[24] Gorham, P. et al. 2002, Nucl. Instrum. 
Meth., A490, 476 

[25] Gorham, P. W., Hebert, C. L., Liewer, 
K. M., Naudet, C. J., Saltzberg, D., & 
Williams, D. 2004, Physical Review Letters, 
93, 041101 

[26] Gorham, P. W., Liewer, K. M., Naudet, 
C. J., Saltzberg, D. P., & Williams, D. 2001, 
AIP Conf. Proc. 579: Radio Detection of 
High Energy Particles, 579, 177 

[27] Gorham, P.W. et al. 2004, in preparation 
for PRD. 

[28] de Gouvela Dal Pino, E.M., Lazarian, A. 
2001, ApJ, 560, 358 

[29] Green, K., Rosner, J.L.D., Suprun, A. & 
Wilkerson, J.F. 2003, Nucl. Instr. Meth. 
A, 498, 256 

[30] Greisen, K. 1956, Prog. Cosmic Ray 
Physics., Suppl. 3. 



[31 

[32; 

[33 
[34 
[35; 

[36 



[37 

[38; 
[39 



[40 

[41 

[42; 
[43; 
[44; 
[45; 

[46; 



Greisen, K. 1966, Phys. Rev. Lett., 16, 748 

Hankins, T.H., Ekers, R.D., O'Sullivan, 
J.D. 1996, MNRAS, 283, 1027 



Haswell, C.A., Tajima, T. 
1992, ApJ, 401, 495 



Sakai, J.-L. 



Hillas, A.M. 1984, Ann. Rev. Astron. As- 
trophys., 22, 425 

Horneffer, A., et al., 2003, in: "Proceedings 
of the 28th ICRC, Tsukuba, Japan", Vol. 
HE 1.5.15, p. 969 

Horneffer, A., et al., 2004, in: "Astronom- 
ical Telescopes and Instrumentation 2004: 
Gravitational Wave and Particle Astro- 
physics detectors", Proc. SPIE, Vol. 5500- 
21 

Hough, J.H., Prescott, J.R. 1970, Nature 
227, 591 

Hough, J.H. 1973, J. Phys. A, 6, 892 

Huege, T., & Falcke, H, Proc. of the 6th 
European VLBI Network Symposium, Ros 
E., Porcas R.W., Lobanov, A. P., & Zensus, 
J.A. (eds.), MPIfR, Bonn, Germany (2002), 
p. 25 

Huege, T., & Falcke, H, Astron. Astro- 
phys., 2003, 412, 19 

Huege, T., & Falcke, H, Astron. Astro- 
phys., 2004, subm. 

Jackson, J.D. 1999, Classical Electrodynam- 
ics, 3rd Edition (Wiley), p. 664. 

Jones, F.C., & Ellison, D.C. 1991, Space 
Sci. Rev., 58, 259 

Kahn, F.D., Lerche, I. 1966, Proc. Roy. Soc. 
London, A-289, 206 

Kamata, K., & Nishimura, J. 1958, Prog. 
Theor. Phys. Suppl. 6, 93 

Learned, J. G., & Pakvasa, S., 1995, As- 
tropart. Phys. 3, 267. 



2G 



[47] Lehtinen, N.G., Gorham, P.W., Jacobson, 
A.R. & Roussel-Dupre, R.A. 2004, Phys. 
Rev., D69, 013008 

[48] Linsley, J. 1986, J. Phys. G, 12, 51 

[49] Litwin, C, Rosncr, R. 2001, Phys. Rev. 
Lett., 86, 4745 

[50] Medvedev, M.V. 2003, Phy. Rev. E, 67, 
045401 

[51] Mucke, A., Protheroe, R.J., Engel, R., 
Rachen, J. P., & Stanev, T. 2003, Astropart. 
Phys., 18, 593 

[52] Nagano, M., & Watson, A. A. 2000, Rev. 
Mod. Phys., 72, 3, 689 

[53] Olhocft, G.R. & Strangway, D.W. 1975, 
Earth Plan. Sci. Lett., 24, 394 

[54] Pelletier, G., Kersale, E. 2000, Astron. As- 
trophys., 361, 788 

[55] Prah, J. H. 1971, M.Phil, thesis, University 
of London 

[56] Protheroe, R.J. 2000, in Topics in cosmic- 
ray astrophysics, ed. M.A. DuVernois, Nova 
Science Publishing: New York, pp 258-298 

[57] Protheroe, R.J., Donea, A.-C, Reimer, A. 
2003, Astropart. Phys., 19, 559 

[58] Protheroe, R.J. & Clay, R.W. 2004, Publ. 
Astron. Soc. Australia, 21, 1 

[59] Protheroe, R.J. 2004, Astropart. Phys. 21, 
415 

[60] Rybicki, G.B., Lightman A.P. 1979, Radia- 
tive Processes in Astrophysics, (John Wiley 
& Sons, New York) 



[61 
[62 



[63 
[64 
[65 

[66 
[67 
[68 
[69 
[70 
[71 

[72 
[73 



Saltzberg, D., Gorham, P., Walz, D. ct 
al. 2001, Phys. Rev. Lett., 86, 2802; 



hep-ex/0011001 



Schieler, H. and the KASCADE and 
LOPES collaborations, 2003, in: Parti- 
cle Astrophysics Instrumentation, Peter W. 
Gorham. (ed.), Proceedings of the SPIE, 
Volume 4858, pp. 41-55 

Sorrell, W.H. 1987, ApJ, 323, 647 

Spencer, R.E. 1969, Nature 222, 460 

Stanev, T., Engel, R., Mucke, A., 
Protheroe, R.J., & Rachen, J.P. 2000, Phys. 
Rev. D, 62, 093005 

Suprun, D., Gorham, P. & Rosner, J. 2003, 
Astropart. Phys., 20, 157 

Takahashi, T. et al., 1994, Phys. Rev. E, 
50, 4041 

Takeda, M., et al. 2003, Astropart. Phys., 
19, 447 

Vcrnov, S.N. et al. 1968, Can. J. Phys. 46, 
241 

Wheeler, J.A., & Feynman, R.P. 1949, 
Rev. Mod. Phys., 21, 425 

Yamamoto, T., Mase, K., Takeda, M., 
Sakaki, N., Teshima, M. 2004, Astropart. 
Phys., 20, 405 

Zas, E., Halzen, F., Stanev, T. 1992, Phys. 
Rev. D, 45, 362 

Zatsepin, G.T. & Kuzmin, V.A. 1966, Sov. 
Phys. JETP Lett., 4, 78 



